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Abstract: 

We present montage, a post-processing nucleosynthesis code that combines a traditional network for 
isotopes lighter than calcium with a rapid algorithm for calculating the s-process nucleosynthesis of 
the heavier isotopes. The separation of those parts of the network where only neutron-capture and 
beta-decay reactions are significant provides a substantial advantage in computational efficiency. We 
present the yields for a complete set of s-process isotopes for a SM©, Z = 0.02 stellar model, as a 
demonstration of the utility of the approach. Future work will include a large grid of models suitable 
for use in calculations of Galactic chemical evolution. 
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1 Introduction 

The asymptotic giant branch (AGB) is the final stage 
of the lives of low- and intermediate-mass stars. Evo- 
lution during the AGB phase is domin ated by thermal 
pulses l|Schwarzschild and Harrnlll965l ). These cyclical 
events consist of brief, intense helium-burning pulses 
interspersed with quiescent hydrogen-burning interpulse 
phases. This complex behaviour produces significant 
nucleosynthesis. Of particular significance for the ori- 
gin of the elements heavier than iron is the s-process. 
Neutrons produced by the "C(a, n)^®0 and ^^Ne(a, n)^ 
reactions are captured on to heavy metal isotopes, in 
particular ^''Fe. Through a sequence of neutron cap- 
tures and beta decays this process can produce iso- 
topes as heavy as lead and bismuth. The AGB and 
s-process nucleosynthesis are well covered in the litera- 
tur e; for ex a mple, see recent reviews of A GB evolution 
by iHerwigj ||2005D and the s-process by iBusso et al.l 
as well as references therein and other contri- 
butions in this issue. 

The recent history of s-process nucleosynthetic cal- 
culations is dominated by the work of Roberto Gallino 
and his collaborators, and is reviewed in detail else- 
where in this volume. Much of this work has been 
motivated by the abundant data provi ded in recent 
years by the analysis of meteoritic grains ()Zinnei]l2008l : 
iLugard [20051 ') and by spectroscopic observations of g i- 
ant stars (e.g. iLambert et alJflQQSl : lAbia et"al ] |2002| ). 
Such calculations are challenging for several reasons. 
One problem is the large computational effort required. 
For stars that undergo thermal pulses the AGB is, by a 
large factor, the most difficult part of the evolutionary 
calculations. In order to reduce their run-time AGB 
structure calculations usually only include the reac- 
tions that generate substantial quantities of energy: 



the pp-chains, CNO cycles and helium-burning reac- 
tions. This allows a small nuclear network to be used. 
For nucleosynthetic calculations, however, a much larger 
network is required. To deal with this problem a post- 
processing technique is commonly employed, where the 
calculation of nuclear burning and mixing is decou- 
pled from stellar structure and evolution calculations. 
The isotopes produced and destroyed via the s-process 
represent the majority of the nuclear network. One 
technique to deal with them is to assume an aver- 
age neutron-absorption cross-section for the heavy iso- 
, topes and, by means of a fictitious partic le, count the 
number of neutrons that they absorb (e.g. lBusso et al.l 
1 19881 1. This then allows those neutrons to be dis- 
tributed within the s-process network to calculate their 
nucleosynthetic effect. Alternatively a restricted s- 
process network can be used, containing a sub-sec tion 
of the relevant isotopes (e.g. iKarakas et al. 20 08|). A 
third approach has been followed by iStraniero et al.l 
(|2006l ). who present an evolutionary model containing 
a full nuclear network c oupled directly to th eir stellar 
evolution code (see also ICristallo et aDl2009l l. 

We describe here a numerical technique that allows 
us to compute the nucleosynthesis of all the relevant 
isotopes at the same time with a reasonable compu- 
tational effort. We split the nuclear network into two 
parts and solve them separately. The lower network 
comprises all elements up to and including potassium, 
and some calcium and scandium isotopes. The upper 
network contains all the heavier isotopes. The code 
for evaluating the lower network is briefly discussed in 
Section [5] and that for the upper network, along with 
the interface between the two, is described in Section|31 
We describe a sample model and give its yields in Sec- 
tion H 
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2 MONSOON and the lower net- 
work 

To model the nucleosynthesis of species less massive 
than calcium we use MONSOON, a post-processing nu- 
cleosynthesis code. It calculates the mixing and burn- 
ing of a user-specified set of isotopes simultaneously 
using a standard relaxation method. As its input MON- 
SOON takes the radius, pressure, temperature, density, 
mixing length and velocity as a function of mass from 
separately calculated stellar structure models. We briefly 
cover the details of monsoon's construc tion here: for 
more information the reader is referred to lCanno 
which describes an early version of the code, and Church 
et al. (in preparation). 

MONSOON comprises a one-dimensional, spherically- 
symmetric model. Material undergoing convection is 
split into two streams, one moving upwards and one 
downwards. In a given timestep material is mixed into 
each up-flowing cell from the cell below, and similarly 
into each down-moving cell from the cell above. A 
gradient in the vertical flow causes material to flow 
between the two cells in the same level. The ratio of 
the cross-sectional areas of the up-flowing and down- 
flowing streams can be adjusted. However, in the model 
computed here we assume that the two streams have 
the same cross-section. 

We adopt a mesh composed of several sections, 
the number depending on the evolutionary phase. For 
main-sequence stars the mesh is placed at geometri- 
cally spaced intervals in mass, increasing from the core 
to envelope. When burning shells develop on the RGB 
and AGB they are treated with two additional sec- 
tions of mesh each, also with geometrically-spaced — 
but closer — meshes. Cores and inter-shell regions are 
treated with meshes spaced at constant intervals in 
mass. As burning shells move out through the star the 
sections of mesh that represent them consume points 
from the region above and release points into the re- 
gion below as necessary. Otherwise the mesh is held 
steady to reduce numerical diffusion. Some special 
treatment is needed to preserve sufficient resolution in 
the inter-shell region on the AGB and this is discussed 
in Section [H] It is important to note that the mesh 
that we adopt does not depend on that used for the 
structural input. This provides a substantial saving in 
computational effort; for example, the structure code 
must place many points in the outer envelope to deal 
with ionisation zones, whereas very little nucleosynthe- 
sis of any interest takes place in the envelope and hence 
very few points are required. Similarly we choose the 
timestep independently of that used to calculate the 
star's structure. 

Nuclear rea ctions are utihsed in the paramete rised 
REACLIB form o f iRauscher and Thielemamil (|2000D . This 
utilises a seven-component fit to reaction rates that 
avoids the need for interpolation in large tables. Mul- 
tiple fits may be made in the case of complex rates, in 
particular those exhibiting resonances. For the models 
presented in this paper we use the recommended rates 
of Illiadis (private communication) . 

To produce a ^^C pocket and hence the required 



source of neu trons in low-mass sta rs we follow the pre- 
scription of l|Karakas et al.|[2008f ). For each thermal 
pulse where third dredge-up occurs a partially-mixed 
zone is artificially inserted at the bottom of the con- 
vective envelope at its point of maximum penetration 
into the star. Instead of trying to model the behaviour 
of the (unknown) process causing this partial mixing, 
we mimic its effect by modifying the hydrogen abun- 
dance to include an exponentially decaying profile just 
below the base of the convective zone. The capture 
of these protons by ^'^C naturally leads to a pocket of 
^^C forming, which burns via the ^^C(Qf,n)^^0 reac- 
tion to produce the dominant neutron source for the 
s-process in low-mass AGB stars. We find that for 
a star of this mass (M ~ SM©) a reasonable quan- 
tity of s-processing is obtained by choosing a zone of 
mass roughly one tenth that of the intershell region 
(,Karakas et al.„2008i '). 



3 Treatment of the s-process 



To tre at the s-process we use the network of lStraniero et aP 
(|2006l ') for isotopes from '''*Ca to ^°^Bi. In this part 
of the network the only reactions that we consider 
are neutron captures, beta decays and electron cap- 
tures. In the conditions found in AGB stars proton- 
and alpha-capture reactions are not signiflcant for these 
isotopes. Decays that occur rapidly are taken to be 
instantaneous. The network is terminated by the re- 
actions 

210-0-/ ,3\207t-)i 

Bi(n,a/3) Pb 



and 



'Po(n,a)^"^Pb 



2i''Po(a)^°«Pb. 



The alpha particles produced by these reactions are 
not fed back into the lower network but are retained 
within the upper network and do not react. This is 
unlikely to lead to significant errors as the s-process 
occurs within helium-rich parts of the star where the 
contribution of the mechanism to the helium abun- 
dance is insignificant. 

In order to increase the speed of the code such that 
it can produce useful results, the burning and mix- 
ing processes in the upper network are decoupled. A 
burning step takes place first, using the neutron abun- 
dances provided by the lower network, and the result- 
ing isotopic abundances are then mixed. Th e mixing 
scheme used is, as in the lower network, that of lCannonI 
( 1993), and again the resulting equations are solved us- 
ing a standard relaxation method. As a result of this 
simplification burning only needs to be calculated at 
a single meshpoint at a time, whereas mixing is only 
calculated for a single isotope at a time. This removes 
the need to invert very large matrices and hence makes 
the problem tractable. It may lead to some small loss 
of accuracy, but only for reactions that occur on a 
similar timescale to the mixing processes. Such reac- 
tions are, in these models, confined to the lower net- 
work, which is why simultaneous mixing and burning 
is employed there. The neutron capture reactions that 
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dominate the upper network are in general much faster 
than convective mixing, whereas most of the f3 decays 
are rather slower. 

For reasons of simplicity the mesh and timestep 
used by the upper network are taken from the lower 
network. Occasionally the upper network may require 
a shorter timestep to converge: in these cases the step 
is broken into a set of several intervals of equal dura- 
tion. 

In order to resolve the ^^C pocket and resulting s- 
process reactions well it is necessary to have good mass 
resolution in the intershell. We find that a good com- 
promise between adequate resolution and code perfor- 
mance can be achieved by inserting a large number of 
points, of the order of 50, when the partially mixed 
zone is included at the point of maximum dredge-up. 
We then allow the code to add and remove points in 
the usual manner, provided that the removal of a point 
does not create a composition discontinuity greater 
than a factor of 1.3 between any of the lower-network 
abundances of the newly adjacent points. This is em- 
pirically found to retain sufficient resolution during the 
interpulse period as the pocket burns radiatively. At 
the peak of the next pulse the intershell convection 
zone reaches the remainder of the ^"^C pocket: the 
convective motions flatten the composition profile and 
hence the points are removed. 

3.1 MONSOON and the upper network 

The interface between MONSOON and the upper net- 
work is effected by means of the species *^Ca and 
*^Ca. These two isotopes are included in both net- 
works. Their molar fractions, h' , in the lower network 
are given, in terms of molar fractions h of species in 
the upper network, by 

/l45Ca - 

where the sum runs over all the species in the upper 
network. In the lower network there are two reactions 
that involve both species: 

**Ca(n,7)*5Ca 

and 

^^Ca(g)*^Ca. 

The fictitious particle g is included to count the num- 
ber of neutron captures that the material in the up- 
per network has undergone. The rate of the second 
reaction is chosen to be very fast, so that any neu- 
tron captured produces a g particle in effect instanta- 
neously. This constrains the abundance of ''^Ca in the 
lower network to be very close to zero at all times, al- 
though it could take a small but finite value in regions 
with a very large neutron exposure. The cross-section, 
(74,4q^ ni of the first reaction is calculated as a weighted 
sum across all the n-capture reactions in the upper net- 
work: 

'^"Ca,n = ^ ^iO-i,n/^44Ca (1) 

where the sum once again runs over all the species in 
the upper network. This allows the number of neutrons 
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Figure 1: The helium burning luminosity of our 
SMq, Zq star is shown as a function of time dur- 
ing the AGB. The time at the start of the plot is 
chosen to precede the onset of oscillations in the 
helium-burning luminosity. The sharp spikes are 
thermal pulses, a total of 45 of which can be seen. 



captured in the upper network to be calculated self- 
consistently, and subsequently used to calculate the 
changes in abundance in the upper network. 

4 Model and yield 

As a sample model we present a SM©, Z = 0.02 
star. The structure model sequence was computed us- 
ing the Monash version of the Mou nt Stromlo Stel- 
lar Structure Programme (monstar, iLattanzi ol 119861 : 
iFrost and Lattanzidll996l ). The opacities used by the 
code have been upd ated. At high temperature s we use 
the OPAL opacities (jlglesias and Roger^ 19961), and at 
low temperatures the opacity tables from Fergus on et al.l 
(2005). We use the Kudritzki and Rcimcrs ( 19781) mass- 
loss law on the first giant branch and during core he- 
liu m burning with t? = 0. 4 , and the AGB mass-loss law 
of IVassfliadis and WoodI (I1993D thereafter. Our ini- 
tial abundances are taken from lAnders and Grevessd 
Cl989). The model was evolved from the pre-main se- 
quence through to the late stages of the AGB where 
numerical instabilities prevented any further evolution. 
The helium-burning luminosity as a function of time 
during the AGB is shown in Figure[T] Dredge-up takes 
place during the leist 35 thermal pulses: the very final 
pulse does not evolve as far as dredge-up and hence is 
not included. At the end of the calculation the model 
has lost a total of 1.77 Mq of material in winds and 
has a hydrogen-exhausted core mass of 0.71 Mq, with 
0.52 Mq remaining in the envelope. Given the mass- 
loss rate of 2.16 x 10~^ Mq yi~^ and interpulse period 
of 3.3 X 10* yr prevailing at the end of the simulation 
this suggests that there are two thermal pulses missing, 
including the final pulse that does not converge. 

We find that we have both a larger final core mass 
a nd more thermal pulses than th e corresponding model 
of lKarakas and Lattanziol (j2007l ) , which was calculated 
with a different version of the same code. This is be- 
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cause we have used updated physical data, including 
a larger mixing length parameter and updated low- 
temperature opacities. These changes lead to a later 
onset of the super- wind, hence a larger core and in- 
creased number of thermal pulses. 

4.1 Nucleosynthesis during a sam- 
ple pulse 

For each pulse with dredge-up we add a partially-mixed 
zone as described in Section [2] This then goes on to 
form a ^'^C pocket through the capture of protons on 
to ^^C. As an example we present the evolution of the 
isotopic a bundances duri n g the twelfth thermal pulse. 
Following iKarakas" et al.l (2008') we add a partially- 
mixed zone of mass 10~^ Mq - chosen to be about 
one tenth of the total intershell mass - at the point of 
maximum dredge-up. Figure [5] shows the effect on the 
most structurally-important isotopes. 

Panel (a) shows conditions 10"^ years after the max- 
imum extent of dredge-up. The partially-mixed zone 
has been added and has partially burnt already to pro- 
duce some ^^N and ^^C. Panel (b) is part- way into the 
interpulse phase. The part of the partially mixed zone 
with the largest abundance of protons has burnt into 
CNO equilibrium, so the most abundant CNO isotope 
is ^"^N. Further in to the star a larger quantity of ^'^C 
has been able to form. Panel (c) is from about 2/3 of 
the way through the interpulse period. The majority 
of the ^''C formed has burnt to ^''O, via the neutron- 
producing alpha-capture reaction. The final panel, (d), 
shows the state of the intershell just before the 13th 
thermal pulse. The ^^C has been almost completely 
burned. 

Figure |3] shows the effect of the ^'^C pocket and 
subsequent neutron production on a selection of s- 
process isotopes. The discontinuity between the abun- 
dances in the envelope and those in the intershell is 
caused by previous s-process episodes, caused by the 
^^^0(0;, n)^'^0 reaction during the interpulse periods and 
the ^^Ne(Q, n)'^^Mg reaction which is eventually acti- 
vated at the pulse peaks. The chemical signatures of 
these previous thermal pulses remain within the in- 
tershell region and are mixed at each pulse peak by 
the partially-overlapping intershell convection zones. 
The isotopes shown were selected to show the effect 
of a moderate neutron exposure on isotopes with a 
range of different formation and destruction pathways 
and neutron-capture cross-sections. ^^Fe is the ini- 
tially most abundant s-process isotope and hence the 
only discernible effect on its abundance is destruction 
by absorption of neutrons. Similarly, ^*Ni is synthe- 
sised in supernovae, not by the s-process: it is also 
only destroyed by neutron absorption. ^^Co is both 
produced and destroyed at different mass co-ordinates 
in the star, depending on the degree of neutron ab- 
sorption. Starting from the majority seed, ^^Fe, ^^Co 
requires only three neutrons to be captured and hence 
is produced in the regions of the pocket with a lower 
neutron exposure. Towards the centre of the pocket 
where the neutron flux is higher it is on average de- 
pleted. The other two isotopes, *^Rb and '^"^Pb, show 
only production, as they lie further up the s-process 



and hence require a large number of neutrons to be 
absorbed by a seed for their production. 

4.2 Comparison with single-network 
approach 

As part of the testing process we created a version of 
MONSOON with an extended network, that went as far 
as ^**Zr. This was used to model the first few thermal 
pulses of the model, as was a version of montage with 
a cut-down s-process network. The adopted trunca- 
tion of the network reduces the runtime for the MON- 
SOON calculation to a manageable level, allowing us to 
compute the models presented here in under a week. 
The nuclear physics input data for the s-process net- 
work in MONTAGE was changed so that the same re- 
action rates were used as in MONSOON. A selected set 
of intershell abundances from the end of the second 
thermal pulse in both codes are shown in Figure O It 
can be seen that, phenomenologically, the abundance 
profiles look very similar but that there are small dif- 
ferences between the abundances calculated, particu- 
larly at the point of highest neutron exposure. Closer 
investigation of the internals of the codes revealed that 
this was due to the different methods used to interpo- 
late reaction rates in temperature: MONSOON evaluates 
the rate for a given temperature from the reaclib for- 
mulae, whereas the s-process network in montage in- 
terpolates in logT within a pre-computed table. This 
leads to a systematically higher reaction rate in mon- 
tage owing to the shape of the rate function at these 
temperatures. We found that the differences in abun- 
dances were consistent with these inconsistencies in 
reaction rates. 



4.3 Yields 

FoUowing lKarakas and Lattanzij l|2007f ) we provide yields 
for the isotopes that we include (Tables [l] and [2]). To 
save space we present only the initial (Solar) value of 
the mass fraction for each isotope, Xi{0), and its aver- 
age value in the ejecta, (Xi). Our model still retained 
some envelope at the point where the evolution termi- 
nated, and the composition of this material is included 
when taking the average. Although we would expect 
two additional thermal pulses to take place whilst the 
remaining envelope was being lost they are unlikely to 
change the yield greatly as they will represent a small 
fraction of the total number of thermal pulses. A more 
sophisticated approach would be to synthesise two ad- 
ditional thermal pulses, as in iKarakas and Lattanziol 
(l2002). 

The yield, commonly defined as 

P. = ^ mt)-Xm^dt, (2) 

can be found aspt = [{Xi) - Xi{0)] AM, where AM = 
2.29 M0 is the total mass lost by the star, including the 
envelope remaining at the end of the calculation. 
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Figure 2: Isotopic abundances for structurally- 
significant isotopes at four points spanning the 
12tli thermal pulse and subsequent interpulse pe- 
riod. Black lines with crosses show hydrogen, red 
lines with plusses ^He, blue lines with stars ^^C, 
magenta lines with squares ^"^C, cyan lines with tri- 
angles ^''N, brown lines with circles and pink 
lines with diamonds ^^Ne. 
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Figure 3: Isotopic abundances for selected s- 
process isotopes at the same times as the corre- 
sponding panels in Figure 2. Black lines with 
crosses show ^^Fe, red lines with plusses ^^Co, blue 
lines with stars ^®Ni, magenta lines with squares 
®^Rb and cyan lines with triangles ^°^Pb. 
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5 Discussion and future work 




Our results demonstrate that we are able to efficiently 
calculate the nucleosynthesis of a large number of s- 
process isotopes, incorporating both the '^^C(a,n)^^0 
and ^^Ne(a, n)^^Mg neutron sources. The nucleosyn- 
thesis calculations presented here took approximately 
one CPU-week of time on a modern desktop computer. 
As can be seen in Figure[2]we considerably over-resolve 
the intershell region. The number of points here could 
most likely be reduced, which would lead to a substan- 
tial saving in runtime. 

Our next target is to extend these results to a large 
grid of models. Ideally Galactic chemical evolution 
calculations need yields at masses separated by about 
0.5 Mq at intermediate masses, possibly lower at the 
lowest masses, and much closer metallicity intervals 
than those currently available (Izzard, private com- 
munication). This will provide a standard, consistent 
set of yields for use in calculation of Galactic chemical 
evolution. 

MONSOON has hitherto always been used with the 
MONSTAR stellar evolution code. D uring testing we at- 
tempt ed to use the evolution code of lStancliffe and Jeffervl 
(|2007l ') to provide output to drive montage. This 
worked very well; a simple approach of writing out ev- 
ery tenth evolution model to the nucleosynthesis code 
input file provided a frequency of data that converged 
well. Within the uncertainties provided by differences 
in reaction rates between the two codes, the light el- 
ement abundances predicted by the two codes agreed 
well. We conclude that montage is a general post- 
processing nucleosynthesis code, which with a small 
amount of effort we should be able to use with any 
stellar evolution code, suitably modified to print out 
the relevant models. This will allow us, for example, 
to quantify the uncertainties in nucleosynthesis pre- 
dictions caused by different evolutionary codes in a 
self-consistent manner. 



Figure 4: Isotopic abundances for selected s- 
process isotopes at the end of the second thermal 
pulse for the codes with a reduced s-process net- 
work. The top panel shows results for MONSOON, 
the bottom panel those for MONTAGE. Black lines 
with crosses show ^^Fe, red lines with plusses ^*Ni, 
blue lines with stars ^^Co, magenta lines with 
squares ^*Rb and cyan lines with triangles ^^Y. 



6 Summary 

The use of a separate network and subtly different 
computational technique for the main nucleosynthe- 
sis and s-process of a thermally-pulsing AGB star has 
a substantial advantage in terms of the computational 
effort required. We have calculated the full s-process 
yields of a 3 Mq star at Solar metallicity, as a demon- 
stration of the utility of the approach. We intend to 
extend this work to a large grid of simulations suitable 
for inclusion in Galactic chemical evolution models, as 
well as by making the code available as a general post- 
processing tool. 
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Table 1: Yields from a 3M0, Z = 0.02 model (see text for details). Continued in Table 2. 



Isotope i (Xi) Xi(Q) Isotope i {Xi) Xi{0) Isotope i (Xi) Xi{0) 



g 


2, 


.999E- 


-04 


2, 


.977E- 


-04 


55 Fe 


6, 


,302E- 


-07 


0, 


■ OOOE-fOO 


94Mo 


7, 


.128E- 


-10 


5 


.829E-10 


p 


6, 


.205E- 


-01 


6, 


.879E- 


-01 


56 Fc 


1, 


.148E- 


-03 


1, 


.167E-03 


95 Mo 


8 


.523E- 


-09 


1 


.014E-09 


3He 


8, 


.963E- 


-05 


8, 


.544E- 


-07 


57 Fc 


3, 


.590E- 


-05 


2, 


.848E-05 


96 Mo 


1 


.260E- 


-08 


1 


.073E-09 


*He 


3, 


.413E- 


-01 


2, 


.931E- 


-01 


58 Fe 


1, 


.095E- 


-05 


3, 


.688E-06 


97Mo 


4, 


.655E- 


-09 


6, 


.210E-10 


7Li 


7, 


.237E- 


-12 


3, 


.550E- 


-20 


60 Fe 


1, 


.090E- 


-07 


0, 


.OOOE+00 


98 Mo 


1, 


.629E- 


-08 


1, 


.585E-09 


12c 


1, 


.743E- 


-02 


3, 


.035E- 


-03 


59 Co 


4, 


.632E- 


-06 


3, 


.444E-06 


ioomo 


1, 


.286E- 


-09 


6 


.456E-10 


13c 


8, 


.063E- 


-05 


3, 


.654E- 


-05 


60 Co 


6, 


.540E- 


-12 


0, 


.OOOE+00 


97Tc 


9 


.579E- 


-13 





.OOOE+OO 


14c 


3, 


.640E- 


-09 


0, 


.OOOE+00 


58 Ni 


4, 


.858E- 


-05 


5, 


.057E-05 


98Tc 


1, 


.053E- 


-13 





.OOOE+OO 


14n 


2, 


.286E- 


-03 


1, 


.106E- 


-03 


59 Ni 


6, 


.309E- 


-07 


0, 


.OOOE+00 


99Tc 


2, 


.032E- 


-09 





.OOOE+OO 


16n 


2, 


.053E- 


-06 


4, 


.367E- 


-06 


60 Ni 


2, 


.092E- 


-05 


2, 


.OOOE-05 


96Ru 


2 


.376E- 


-10 


2, 


.574E-10 


16o 


8, 


.960E- 


-03 


9, 


.602E- 


-03 


61 Ni 


1, 


.493E- 


-06 


8, 


.804E-07 


98Ru 


8, 


.579E- 


-11 


8, 


.948E-11 


i^O 


1, 


.906E- 


-05 


3, 


.891E- 


-06 


62 Ni 


3, 


.731E- 


-06 


2, 


.843E-06 


99Ru 


6 


.965E- 


-10 


6 


.106E-10 


18o 


1, 


.520E- 


-05 


2, 


.169E- 


-05 


63 Ni 


9, 


.193E- 


-08 


0, 


.OOOE+00 


iooRu 


7, 


.636E- 


-09 


6 


.119E-10 


19p 


5, 


.940E- 


-06 


4, 


.055E- 


-07 


64Ni 


1, 


.662E- 


-06 


7, 


.438E-07 


101 Ru 


2, 


.346E- 


-09 


8, 


.339E-10 


20Ne 


1, 


.648E- 


-03 


1, 


.621E- 


-03 


63 Cu 


9, 


.166E- 


-07 


6, 


.103E-07 


102 Ru 


9 


.467E- 


-09 


1, 


.565E-09 


21 Ne 


5, 


.878E- 


-06 


4, 


.131E- 


-06 


65 Cu 


6, 


.137E- 


-07 


2, 


,806E-07 


104 Ru 


1 


.672E- 


-09 


9 


.446E-10 


22Ne 


3, 


.954E- 


-03 


1, 


.303E- 


-04 


64 Zn 


1, 


.154E- 


-06 


1, 


.046E-06 


103 Rh 


2 


.630E- 


-09 


9 


.315E-10 


23Na 


1, 


.334E- 


-04 


3, 


.342E- 


-05 


66 Zn 


1, 


.007E- 


-06 


6, 


.186E-07 


104pd 


5 


.413E- 


-09 


4, 


.211E-10 


24 Mg 


5, 


.184E- 


-04 


5, 


.153E- 


-04 


67 Zn 


1, 


.787E- 


-07 


9, 


.229E-08 


105 Pd 


2 


.176E- 


-09 


8 


.522E-10 


28 Mg 


1, 


.070E- 


-04 


6, 


.773E- 


-05 


68 Zn 


8, 


.792E- 


-07 


4, 


.283E-07 


106pd 


7 


.181E- 


-09 


1 


.053E-09 


26 Mg 


1, 


.349E- 


-04 


7, 


.768E- 


-05 


™Zn 


1, 


.448E- 


-08 


1, 


.458E-08 


107pd 


1 


.223E- 


-09 





.OOOE+OO 


27A1 


6, 


.327E- 


-05 


5, 


.804E- 


-05 


69 Ga 


1, 


.118E- 


-07 


4, 


.019E-08 


108 Pd 


8 


.646E- 


-09 


1 


.038E-09 


2«Si 


6, 


.495E- 


-04 


6, 


.537E- 


-04 


71 Ga 


9, 


.499E- 


-08 


2, 


,744E-08 


iiopd 


1 


.074E- 


-09 


4, 


.515E-10 


2»Si 


3, 


.690E- 


-05 


3, 


.429E- 


-05 


70 Ge 


1, 


.676E- 


-07 


4, 


.397E-08 


107 Ag 


6 


.344E- 


-10 


6 


.925E-10 


30Si 


2, 


.655E- 


-05 


2, 


.354E- 


-05 


72 Ge 


2, 


.059E- 


-07 


6, 


.045E-08 


109 Ag 


2 


.584E- 


-09 


6 


.554E-10 


31p 


1, 


.112E- 


-05 


8, 


.163E- 


-06 


73 Ge 


5, 


.901E- 


-08 


1, 


.745E-08 


108Cd 


5 


.406E- 


-11 


3 


.974E-11 


32g 


3, 


.929E- 


-04 


3, 


.962E- 


-04 


74Ge 


2, 


.678E- 


-07 


8, 


.276E-08 


iiocd 


6 


.640E- 


-09 


5 


.680E-10 


33g 


3, 


.840E^ 


-06 


3, 


.225E- 


-06 


76 Ge 


1, 


.807E- 


-08 


1, 


.816E-08 


iiiCd 


2 


.701E- 


-09 


5 


.874E-10 


34g 


1, 


.913E- 


-05 


1, 


.868E- 


-05 


75 As 


2, 


.930E- 


-08 


1, 


.263E-08 


112Cd 


8, 


.943E- 


-09 


1, 


.117E-09 


36S 


1, 


.499E- 


-07 


9, 


.389E- 


-08 


76 Se 


7. 


.521E- 


-08 


1, 


.279E-08 


113Cd 


2, 


.829E- 


-09 


5, 


.709E-10 


36 CI 


2, 


.501E- 


-06 


2, 


.535E- 


-06 


77Se 


3, 


.546E- 


-08 


1, 


.094E-08 


114Cd 


1, 


.231E- 


-08 


1, 


.354E-09 


36 CI 


5, 


.508E- 


-09 


0, 


.OOOE-fOO 


78 Se 


1, 


.268E- 


-07 


3, 


.442E-08 


iioCd 


1, 


.461E- 


-09 


3 


.592E-10 


37C1 


1, 


.792E- 


-06 


8, 


.554E- 


-07 


79 Se 


2, 


.608E- 


-08 


0, 


.OOOE+00 


ll3ln 


2 


.108E- 


-11 


2, 


.306E-11 


36 Ar 


7, 


.620E- 


-05 


7, 


.748E- 


-05 


80Se 


2, 


.554E- 


-07 


7, 


.525E-08 


iisin 


2 


.788E- 


-09 


5, 


.224E-10 


38 Ar 


1, 


.580E- 


-05 


1, 


.540E- 


-05 


82 Se 


1, 


.331E- 


-08 


1, 


.342E-08 


li4Sn 


6 


.995E- 


-11 


7, 


.450E-11 


39 Ar 


1, 


.962E- 


-10 


0, 


.OOOE+OO 


79Br 


1, 


.790E- 


-08 


1, 


.243E-08 


ilSSn 


3 


.643E- 


-11 


3 


.866E-11 


40 Ar 


6, 


.522E- 


-08 


2, 


.535E- 


-08 


81 Br 


3, 


.695E- 


-08 


1, 


.239E-08 


iiesn 


1, 


.530E- 


-08 


1, 


.672E-09 


39k 


3, 


.528E- 


-06 


3, 


.473E- 


-06 


80 Kr 


1, 


.166E- 


-08 


2, 


.605E-09 


ll7Sn 


5, 


.088E- 


-09 


8, 


.902E-10 


40k 


4, 


.OlOE- 


-08 


5, 


.242E- 


-09 


81 Kr 


4, 


.839E- 


-09 


0, 


. OOOE+00 


llSSn 


2, 


.116E- 


-08 


2, 


.834E-09 


4lK 


4, 


.194E- 


-07 


2, 


.634E- 


-07 


82 Kr 


1, 


.222E- 


-07 


1, 


.377E-08 


ll9Sn 


6, 


.378E- 


-09 


1 


.013E-09 


40Ca 


5, 


.914E- 


-05 


5, 


.996E- 


-05 


83 Kr 


5, 


.153E- 


-08 


1, 


.396E-08 


l20Sn 


3 


.009E- 


-08 


3 


.879E-09 


«Ca 


2, 


.388E- 


-08 


0, 


.000E4 


-00 


84Kr 


2, 


.983E- 


-07 


7, 


.035E-08 


i22Sn 


1 


.634E- 


-09 


5 


.604E-10 


42Ca 


5, 


.076E- 


-07 


4, 


.200E- 


-07 


85 Kr 


1, 


.039E- 


-10 


0, 


. OOOE+00 


i24Sn 


7 


.028E- 


-10 


7 


.117E-10 


«Ca 


1, 


.065E- 


-07 


8, 


.982E- 


-08 


86 Kr 


6, 


.699E- 


-08 


2, 


.197E-08 


121Sb 


2 


.625E- 


-09 


5 


.343E-10 


44Ca 


1, 


.532E- 


-06 


1, 


.476E- 


-06 


85 Rb 


5, 


.745E- 


-08 


1, 


.133E-08 


123Sb 


7. 


.322E- 


-10 


4 


.038E-10 


■leca 


5, 


.266E- 


-09 


2, 


.959E- 


-09 


87Rb 


2, 


.857E- 


-08 


4, 


.471E-09 


122 Tc 


3 


.312E- 


-09 


3 


.915E-10 


48Ca 


1, 


.439E- 


-07 


1, 


.443E- 


-07 


86 Sr 


9, 


.243E- 


-08 


5, 


.098E-09 


i23Te 


1 


.138E- 


-09 


1 


.367E-10 


«Sc 


5, 


.754E- 


-08 


4, 


.069E- 


-08 


87Sr 


6, 


.321E- 


-08 


3, 


,605E-09 


124Tc 


6 


.625E- 


-09 


7 


.370E-10 


«Ti 


2, 


.492E- 


-07 


2, 


.302E- 


-07 


88 Sr 


9, 


.072E- 


-07 


4, 


,342E-08 


125 Tc 


3 


.028E- 


-09 


1 


.108E-09 


47TJ 


2, 


.142E- 


-07 


2, 


.146E- 


-07 


90 Sr 


2, 


.532E- 


-10 


0, 


■ OOOE+OO 


126Tc 


1 


.373E- 


-08 


2 


.972E-09 


«Ti 


2, 


.156E- 


-06 


2, 


.216E- 


-06 


89 Y 


1, 


.827E- 


-07 


1, 


.086E-08 


128Tc 


6 


.148E- 


-09 


5 


.071E-09 


49Ti 


2, 


.356E- 


-07 


1, 


.686E- 


-07 


90 Y 


6, 


.435E- 


-14 


0, 


■ OOOE+OO 


130Tc 


5 


.387E- 


-09 


5 


.514E-09 


60Ti 


2, 


.625E- 


-07 


1, 


.689E- 


-07 


90 Zr 


1, 


.589E- 


-07 


1, 


.355E-08 


127j 


4 


.176E- 


-09 


2 


.953E-09 


SOy 


8, 


.952E- 


-10 


9, 


.402E- 


-10 


91 Zr 


4, 


.272E- 


-08 


2, 


.988E-09 


129j 


2, 


.874E- 


-11 





.OOOE+OO 


Sly 


3, 


.845E- 


-07 


3, 


.826E- 


-07 


92 Zr 


6, 


.210E- 


-08 


4, 


.617E-09 


128Xc 


3 


.588E- 


-09 


2 


.979E-10 


60cr 


7, 


.274E- 


-07 


7, 


.643E- 


-07 


93 Zr 


1, 


.867E- 


-08 


0, 


.OOOE+OO 


129XC 


4, 


.828E- 


-09 


3 


.748E-09 


8lCr 


6, 


.988E- 


-10 


0, 


.OOOE-fOO 


94 Zr 


6. 


.944E- 


-08 


4, 


.781E-09 


130XC 


7. 


.566E- 


-09 


6 


.OlOE-10 


S2Cr 


1, 


.514E- 


-05 


1, 


.533E- 


-05 


96 Zr 


4, 


.960E- 


-09 


7, 


.866E-10 


i3iXe 


5 


.410E- 


-09 


3 


.020E-09 


63 cr 


1, 


.753E- 


-06 


1, 


.771E- 


-06 


93Nb 


1, 


.618E- 


-09 


1, 


.678E-09 


i32Xe 


1, 


.622E- 


-08 


3 


.718E-09 


64cr 


6, 


.234E- 


-07 


4, 


.492E- 


-07 


94 Nb 


3, 


.326E- 


-13 


0, 


.OOOE+OO 


i34Xe 


2 


.066E- 


-09 


1, 


.390E-09 




1, 


.324E- 


-05 


1, 


.339E- 


-05 


92 Mo 


8, 


.647E- 


-10 


9, 


.152E-10 


l36Xe 


1, 


.146E- 


-09 


1, 


.148E-09 




6, 


.855E- 


-05 


7, 


.114E- 


-05 


93 Mo 


5, 


.236E- 


-12 


0, 


.OOOE+OO 
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Table 2: Yields from a SMq, Z = 0.02 model continued from Table 1. 



Isotope i (Xj) Xi(0) Isotope i (Xj) Xi{0) Isotope i {Xj} Xi(0) 



133 Cs 


3, 


.064E- 


-09 


1, 


.289E- 


09 


157 Gd 


2, 


.898E- 


10 


2, 


.172E- 


-10 


184 -y^ 


7, 


.342E- 


-10 


1 


.987E- 


10 


135 Cs 


9, 


.940E- 


■10 


0, 


.000E4- 


00 


158 Gd 


6, 


.980E- 


10 


3, 


.470E- 


-10 


186 -\y 


4 


.370E- 


-10 


1 


.873E- 


-10 


137 Cs 


4, 


.150E- 


■11 


0, 


.000E4- 


-00 


160 Gd 


3, 


.439E- 


10 


3, 


.093E- 


10 


185 Re 


2 


.113E- 


-10 


9 


.477E- 


-11 


l34Ba 


4, 


.353E- 


■09 


3, 


.862E- 


10 


159 Tb 


3, 


.224E- 


10 


2, 


.558E- 


10 


187Re 


1 


.881E- 


-10 


1 


.603E- 


-10 


135 Ba 


2, 


.333E- 


09 


1, 


.061E- 


-09 


160 


1, 


.594E- 


-10 


3, 


.980E- 


11 


186 Os 


2 


.085E- 


-10 


5 


.186E- 


11 


136 Ba 


1, 


.229E- 


08 


1, 


.274E- 


09 


i6lDy 


3, 


.616E- 


-10 


3, 


.236E- 


10 


1870s 


1 


.063E- 


-10 


5 


.280E- 


-11 


l37Ba 


1, 


,045E- 


-08 


1, 


.835E- 


09 


i62Dy 


6, 


.797E- 


-10 


4, 


.393E- 


10 


188 Os 


6 


.752E- 


-10 


4, 


.412E- 


-10 


138 Ba 


7, 


.528E- 


08 


1, 


.180E- 


08 


163 Dy 


4, 


.174E- 


10 


4, 


.315E- 


10 


189 Os 


5 


.666E- 


-10 


5 


.370E- 


-10 


139 La 


8, 


.706E- 


-09 


1, 


.656E- 


-09 


164 Dy 


9, 


.496E- 


10 


4, 


.915E- 


10 


190 Os 


1 


.191E- 


-09 


8 


.851E- 


10 


140 Ce 


1, 


.962E- 


-08 


3, 


.797E- 


09 


163 jjo 


7, 


.238E- 


11 


0, 


.000E+ 


00 


192 0s 


1 


.385E- 


-09 


1 


.389E- 


-09 


i«Ce 


7. 


.390E- 


-10 


4, 


.823E-- 


10 


165 Ho 


4, 


.659E- 


10 


3, 


.838E- 


10 


194 0s 


1 


.072E- 


-15 





.OOOE+00 


141 Pr 


2, 


.081E- 


09 


6, 


.393E- 


10 


164 Er 


5, 


.050E- 


11 


1, 


,771E- 


11 


191 Ir 


1 


.181E- 


-09 


1 


.201E- 


09 


142 Nd 


4, 


.341E- 


-09 


8, 


.555E- 


10 


i66Er 


5, 


.211E- 


10 


3, 


.741E- 


10 


193 Ir 


1 


.979E- 


-09 


2 


.039E- 


09 


143 Nd 


9, 


.027E- 


-10 


3, 


.868E- 


-10 


l67Er 


3, 


.257E- 


-10 


2, 


.571E- 


-10 


192 Ft 


1, 


.630E- 


-10 


5, 


.273E- 


-11 


144 Nd 


2, 


.363E- 


-09 


7, 


.611E- 


■10 


168 Er 


6, 


.655E- 


-10 


3, 


.020E- 


-10 


193 Pt 


1, 


.465E- 


-12 





.OOOE+00 


145Nd 


5. 


.489E- 


■10 


2, 


.673E- 


-10 


170 Er 


3, 


.177E- 


-10 


1, 


.699E- 


-10 


194pt 


2 


.584E- 


-09 


2, 


.247E- 


■09 


146Nd 


2, 


.006E- 


-09 


5, 


.573E- 


-10 


169 Tm 


2, 


.712E- 


-10 


1, 


.716E- 


-10 


195pt 


2 


.334E- 


-09 


2, 


.321E- 


■09 


148Nd 


3, 


.658E- 


-10 


1, 


.893E- 


-10 


i70Yb 


1, 


.397E- 


-10 


3, 


.399E- 


11 


196 Ft 


2 


.467E- 


-09 


1, 


.739E- 


-09 


ISONd 


1, 


.773E- 


-10 


1, 


.879E- 


-10 


171 Yb 


3, 


.077E- 


-10 


1, 


.603E- 


-10 


188 Ft 


5 


.222E- 


-10 


5, 


.013E- 


-10 


I46pm 


2, 


.278E- 


■15 


0, 


.OOOE+00 


172 Yb 


5, 


.841E- 


-10 


2, 


.469E- 


-10 


197 Au 


1, 


.163E- 


-09 


1, 


.003E- 


-09 


144 Sm 


2, 


.883E- 


-11 


3, 


.063E- 


-11 


173 Yb 


3, 


.432E- 


-10 


1, 


.828E- 


-10 


198Hg 


8, 


.556E- 


-10 


1, 


.914E- 


-10 


146 Sm 


2, 


.254E- 


■13 


0, 


.OOOE+00 


174Yb 


1, 


.165E- 


-09 


3, 


.627E- 


-10 


199Hg 


6 


.248E- 


-10 


3 


.256E- 


■10 


l47Sm 


3, 


.150E- 


■10 


1, 


.513E- 


-10 


176Yb 


2, 


.441E- 


-10 


1, 


.465E- 


-10 


200Hg 


1, 


.413E- 


-09 


4 


.480E- 


-10 


148 Sm 


4, 


.837E- 


■10 


1, 


.148E- 


-10 


175 Lu 


2, 


.720E- 


-10 


1, 


.653E- 


-10 


20lHg 


6 


.436E- 


-10 


2, 


.569E- 


-10 


149 Sm 


2, 


.002E- 


■10 


1, 


.411E- 


-10 


176 Lu 


2, 
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